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ABSTRACT 

The heat transfer and pressure drop characteristics of water flow- 
ing in a vertical annulus under counter flow conditions were investigat- 
ed* The heat exchanger used consisted of a straight tube in tube test 
section with water as the heat transfer medium in both tubes* Convec- 
tion coefficients of heat transfer between the inner tube and the fluid 
flowing in the annulus were determined from the measured inner tube wall 
and average fluid bulk temperature and the energy transferred to the 
water flowing in the annulus* 

The friction factor data indicate that the proper choice of dia- 
meter to be employed in evaluating the Reynolds number is the equi- 
valent diameter, de * C show that the transition from 

laminar to turbulent flow occurs at Re ^ 2000 as it does for flow in- 
side tubes. The results for the laminar region are in excellent agree- 
ment with the theoretical predictions. Data in the turbulent region 
correlate well with that in the literature for commercial pipes. 

The heat transfer data cover a range of Reynolds numbers from 200 
to 5000. Two different types of runs were madej equal capacitance, 

(equal mass rate of flow in both tubes for this particular investiga- 
tion) and constant wall temperature. The data indicate that the bulk 
temperature measurements were not reliable. However, since the measured 
temperature differences were reliable, it was possible to correlate both 
laminar and turbulent flow data by one empirical equation. The bulk tem- 
perature distribution for the constant wall temperature run*^ along the 
annulus length was found to have a large effect on the heat transfer 
coefficients* It is not linear and should be taken into account. The 
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data were best correlated by evaluating fluid properties at an average 
film temperature rather than an average bulk temperature as recommend- 
ed in most of the literature. 
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I Introduction. 



The problem of heat transfer in annular spaces has not been ade- 
quately solved and is complicated by the many forms in which the em- 
pirical correlation equations are presented, f 3^ , ^4^ , 

Although a number of equations relating convection heat transfer co- 
efficients to system parameters have been presented in recent years, 
[i] . [2] , these provide correlations only for selected ranges 

of the pertinent variables. The development of proper relations for 
predicting annular heat transfer is made difficult, not only by an 
almost complete lack of data in the laminar region » *^ut also by 

insufficient data covering both laminar and turbulent flow in the same 
annulus. 

It has been the objective of this investigation to obtain (1) fric- 
tion data for water flowing in the annulus of a vertical double tube 
heat exchanger over a range of Reynolds number 180 to 8000 and (2) heat 
transfer data for the same annulus over a range of Reynolds number 200 
to 5000 which was the limit of .the equipment. The purpose of the fri- 
ction data was to determine the proper characteristic length to be 
utilized in evaluating the Reynolds number. The heat transfer results 
were treated in a manner to bring out their comparison with existing 
empirical correlations, and to establish a correlation which would cover 
the transition Reynolds number range. 



(*) Numbers in brackets refer to the Bibliography on page 19 . 
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II. Description of Equipment. 

Heat Exchanger . The heat exchanger proper consisted of a single 
vertical 1.125 inch 0. D. x 0.995 inch I. D. 99.9% copper tube inside 
a 1.372 inch 0. D. x 1.242 inch I. D. 99.9% copper tube 142 inches 
long. The annular diametric clearance was 0.117 inches and was main- 
tained by a helical arrangement of 36 stainless steel ball bearingSj, 

0.093 inch diameter, imbedded in the inner tube wall and fastened by 
solder. Details of construction are shown in Figures 1 and essential 
dimensions in Table 1 

Certain aspects are worthy of comment. The annular space was made 
sufficiently small so that reasonable annulus water velocities might be 
obtained at low mass rates of flow. Inlet and outlet header spaces were 
made sufficiently large to permit proper measurement of the bulk tempera- 
tures. Copper-cons tantan immersion thermocouples were used to measure 
the inlet and outlet temperatures of both streams. A switch was built 
which enabled both absolute temperatures and temperature differences to 
be read. A schematic wiring diagram of the switching arrangement is 
shown in Figure 2. The circuit was designed to eliminate the interference 
from the circulating circuits existing between all the thermocouples. 

On the outer annulus tube 36 copper-constantan thermocouples were anchor- 
ed by soldering in 1/32 inch holes drilled half way through the tube 
wall. These were placed in groups of three, equally spaced around a 
diameter, at 11 inch intervals down the length of the tube. The in- 
side wall temperature of the inner tube was measured by means of a 
^traveling** thermocouple. Several designs were used and Figures 3 and 
'4 show the details of construction. A Victory Corporation 2000 ohm 
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glass enclosed thermistor was used for one design of the "traveling" 
thermocouple. It was calibrated in water with all leads in place and 
the calibration curve is shown in Figure 15, 

Pressure taps were drilled in the outer tube at a distance of 21 
inches from each end giving a distance of 100 inches between taps. Both 
mercury and water manometers were used to measure pressure drop. For 
reading the small pressure differentials a Fisher cathetometer was used 
which read to the nearest 0.01 centimeter. 

Auxiliary Equipment . The overall test assembly is shown in Figure 
5, A large hot water storage tank provided a heat reservoir of uniform 
temperature. The water was not changed during the investigation, there- 
fore assuring a minimum of entrapped gasses and particles in the circuit 
Fischer and Porter Flowrator meters were used to measure the flow rates 
for both hot and cold systems. These also served to indicate stability 
of the flow. The 36 outer wall thermocouple readings were taken with a 
Minneapolis - Honeywell 48 position recording potentiometer. These read 
ings were recorded to a precision of + 0.2°F. All other thermocouples 
were read with a portable precision hand balance Rubicon potentiometer 
which was accurate to the nearest +0,1 ®F. The thermistor resistance 
was read with a wheatstone bridge balancing arrangement and was accurate 
to the nearest +0.1 ’F. 

The hedt exchanger proper was wrapped with 2 inch width asbestos 
cloth and then covered with 1 1/2 inch thick pipe lagging. 
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IIIo Experimental Procedure. 



Friction : To obtain isothermal friction data, pressure drop 

measurements were made of the fluid flowing in the annular space at 
various flow rates. The pressure taps were so located to give an entrance 
length of 180 equivalent diameters. Pressure differences of less than 
50 inches of water were determined with a water manometer using a catheto* 
meter to secure accurate readings. For pressure differences of 50 inches 
of wat:er or more a mercury manometer was used. The water temperature 
for most of the runs was essentially room temperature. Five check runs 
were made at a temperature of 100®F. For low flow rates an electric 
timer and a weighing tank weVe used to determine the mass rate of flow. 
These measurements were used to check the calibration of the flow meters* 
At low flow rates the flow meter measurements varied as much 6% with 
the flow rate obtained by weighing, but the overall variation was less 
than + 2%. 

Heat Transfer : The procedure for the equal capacitance luns and the 

constant wall temperature runs was essentially the same. For equal capa- 
citance, the mass flow rates in the annulus and inner tube were kept the 
same since the specific heat for water was constant over the temperature 
range employed. Flow rates for the annulus were varied to cover a 
Reynolds number range of 200 to 5000. 

For constant wall temperature runs the inner tube mass flow rate 
was held constant at as high a value as permitted by the apparatus, 
corresponding to a Reynolds number of about 11,000 for all runs. For 
the annulus, the mass flow rate was varied to correspond to a range of 
Reynolds number of 200 to 2800. 
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The hot water temperature was controlled by rheostats on the 
heaters and the energy input measured by a voltmeter and ammeter. 

For the hot water circulating system sufficient cold make-up water was 
added to keep a constant water level in the reservoir. 

Data for the runs were taken after five to ten minutes of steady 
state operation of the system. Steady state conditions were denoted by 
the constancy of flow rates and temperatures and were usually reached 
within an hour after changes in operating conditions had been made. Re- 
peat data recordings were made of all temperatures except those of the 
"traveling" thermocouple. 

The inner tube wall temperatures were determined by moving the 
"traveling" thermocouple up or down the tube in 12 inch increments. 
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IV. Experimental Results. 

The friction results are shown in Figure 6 and tabulated in Table 
2. The table also presents the temperatures, pressure differences and 
flow rates for the individual runs. 

The heat transfer results are shown in Figures 7, 8, 9, 10 and 17 
and tabulated in Table 3. The table also presents the temperatures, flow 
rates and all other data used to calculate the system parameters as well 
as the parameters themselves. Inner wall temperatures were not correct- 
ed for thermocouple error as calculations of the temperature difference 
introduced by the exposed wires, using Schneider's flat plate analogy 
C63, showed this effect to* be insignificant. 

The inner wall temperatures were corrected according to laminar 
or turbulent flow in the inner tube for the physical distance of the 
thermocouple from the inner wall. This distance was 0.02 inches and 
amounted to 4% of the tube radius. Figure 14 shows Martinelli's [163 
and Eckert and Drake's [153 temperature profiles for turbulent and 
laminar flow respectively. A value of 0.625 was used for turbulent 
flow calculations and a value of 0.15 for laminar flow calculations. 

A correction to account for the temperature drop across the inner tube 
wall was also made. The value of thermal conductivity, k, for pure 
copper of ,220 BTU/hr. ft. ‘“F was used. This value is probably slightly 
high as the tube is only 99.9% pure copper. 

The average bulk temperatures of the annular fluid for constant 
wall temperature runs were corrected to account for the temperature 
distribution with annulus length. The three thermocouple readings at 
each of the 12 outside annulus wall positions were averaged and a 
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dimensionless temperature, t^ - tac /^ca » was computed for 

each position. Simpsons Rule was used to calculate an average bulk 
temperature for each run. The average inner annulus wall tempera- 
tures were determined by an arithmetic average of the end wall tempera- 
tures. 

For the equal capacitance runs both the inner and outer annulus wall 
temperature distributions were considered linear and the average wall 
and average bulk temperatures determined on that basis. 

The system parameters are based on a film value temperature^ t^ ^ 
defined as 

U - t',n + tbc»v fa (1) 

The data were also correlated using bulk temperatures and an aver- 
age bulk temperature as a basis for calculating system parameters. 
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V, Sources of Error. 



The temperature measurement technique offered the prfncijal chance 
for error and all numerical results were a direct consequence of these. 

To minimize these errors, all thermocouples for bulk temperature measure- 
ments were similar and all manufactured thermocouples came from the same 
spool of wire and were made by the same method. A circuit arrangement, 
Figure 2, was built to eliminate circulating current effects. Measured 
differences rather than absolute temperature measurements were used 
whenever possible. All readings were made on the same potentiometer 
and referenced against the same standard cell. 




As shown in Figure A, the sum of AtTii + should equal the sum 

^^34 • These measured temperature differences were 

totaled and found to vary an average of 2 5%, Tils indicates the circuit 
design eliminated the circulation errors. 

The sensitivity of the potentiometer was + .001 millivolt or 
approximately + 0.04®F. Conversion of millivolts to degrees Fahren- 
heit was done from Leeds and Northrup Conversion Tables for Thermocouples. 
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Flow rate variation was a second potential source of error. The 
flow meter readings were checked during the friction test by weighing 
the mass of water for a given time. The flow meters are considered 
accurate to + 2%. At least five readings were taken for each run and 
they did not vary more than 1%. 

Further sources of error which were determined to be insignificant 
included (a) expansion or contraction of the annulus tubes with tempera 
ture fluctuations, (b) variation of the temperature during the time re- 
quired for observation, (c) heat loss along the thermocouple wire and 
(d) the bi-metallic thermoelectric effect at the connector terminals. - 
The latter connectings were designed and located to insure that both 
connections for each thermocouple lead would remain at the same tempera 
ture thus assuring no effect on the results from this extra junction. 

An uncertainty analysis of the above effects on the system para- 
meters indicates that the repeatability of experimental data should 
be within + 10%. 



9 



VI, Discussion of Results. 



Friction; The isothermal friction data shown in Figure 6 are plot- 
ted as friction factor versus Reynolds number, based on an equivalent 
diameter. The change from laminar to turbulent flow is seen to occur 
approximately at a Reynolds number of 2000 which corresponds to the 
critical value for transition of flow inside tubes Do] . As it is 
desirable for comparison purposes to have the critical value of Reynolds 
number which conforms with that of flow in tubes, it was decided the 
equivalent diameter was the preferred characteristic defining length. 

As shown by Lamb [7] , besides the usual relations for friction 
factor. 



and Reynolds number. 






( 2 ) 



Re = A& - 

il 7 



(3) 



for flow of fluids in an annulus an additional geometry parameter, 



0 (d./J.) = 



(i- d./j>f 

I + 






( 4 ) 



is needed 



The limits of being 

wK«vv e/./c/v-O 







'^1 — vll- 



§ ♦ • 




as for flow in tubes j and 



as for flow between parallel plates. The functional relationship 
between the parameters for laminar flow is theoretically found to be 

(5) 

For the annulus used in this investigation 

diChu^) = 

giving a theoretical prediction of the friction factor for laminar 
flow as 

'It * 

From F Lgure 6 it is seen that the experimental values of the friction 
factor are within + 5% of the theoretical values over a Reynolds 
number range of 200 to 2000, 

In the turbulent flow region^ (Re > 2000), the data were plotted 
and correlated by the equations 

( 7 ) 

These data plot 12% higher than the proposed equation of Knudsen and 
Katz C8D of; 

( 8 ) 

and 6«5% higher than the proposed Davis CO relationship ; 
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(9) 



The correct function of the Reynolds number to use appears to 
be (Re) as the slope of the line in Figure 6 is -0.2. Equation 



(8) expresses the friction factor by taking the average of a number 
of investigators* experimental results {^8^ and neglects to account for 
any function of diameter. Knudsen and Katz C.83 predict a 35% varia- 
tion from equation (8), depending on the ratio of ). Only 

for larger (Jv/Jj ) ratios, above 3, may equation (8) be used with 
small error. It is, therefore, recommended that the Davis equation be 
used, especially for small values of ). It is noted that the 

turbulent data fall very close to the line for flow Inside commercial 
pipes. 

Heat Transfer . Figure 7 is a plot of all the laminar flow data 
correlated by using the Chen, Hawkins and Solberg t^l^ equation based 
on the arithmetic mean bulk temperature of the fluid and the viscosity, 
, evaluated at the inner annulus wall temperature t,^ t 



was used and a plot is presented in Figure 8. The turbulent data fall 
50% below the empirical curve. 

The lack of correlation between these data and the empirical equations 
is obvious. The main explanation lies in the fact that the measured bulk 




The data fall an average of 30% below the empirical curve. 

For the turbulent equal capacitance runs. Monrad and Felton's [12J 



equation based on the average fluid bulk temperatures. 




( 11 ) 
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temperatures were not true temperatures. The mixing chambers did not 
perform properly and further, the inner pipe conducted heat into both 
annuli mixing chambers, which would also contribute to bulk temperature 
error. This is substantiated by comparing the hot and cold bulk tempera- 
ture readings at the ends of the annulus. The cold water inlet tempera- 
ture which should remain fairly constant rises and falls with a hot 
water temperature rise and fall. 

The hot water bulk temperature measurements are considered better 
than the cold water bulk temperature measurements because the thermo- 
couples could be placed nearer the center of the tube. The cold water 
bulk thermocouples could not be placed closer than the outside radius 
of the inner tube and, therefore, could be considerably influenced by 
conduction from the hot tube as well as not actually "seeing" the in- 
let or outlet cold water temperatures. The inner and outer annulus wall 
temperatures are considered good as well as the bulk temperature differ- 
ences. 

Another possible reason for the discrepancies with these correla- 
tions is due to the effect of natural comection. There are no corre- 
lations in the literature that cover the flow conditions for the particu- 
lar annulus investigated. In general, natural convection is expected to 
influence the heat transfer behavior in an amount of + 15% depending on 
its relation to the external forces Cl3] . For the Chen^ Hawkins and 
Solberg correlations with the Grashof number raised only to the 0o05 power^ 
practically all of the free convection effects are damped onto In the 
Monrad and Felton correlation no free convection effects are included 
since^ with Reynolds numbers in excess of 10,000, the inertia forces 
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completely override the buoyant forces. 

It was found that all of the data of this investigation, in- 
cluding the transition Reynolds number range, could be correlated by 
the equation 

'iy X't5' .OtpS' ■S’lle 

Nu = 0.093 (Re) (Pr) (Gr) (^c/t) (V«/,) (12) 

a plot of which is presented in Figure 10. This equation correlated 
all data with an average variation of + 107o and a maximum variation of 
207o. The Sieder and Tate (ll73 viscosity correlation, (<*^/ui), is not 
included as it is recognized that the present data do not cover a suf- 
ficient variation in temperature gradient to test this relationship, 
and also previous investigations have not been sufficiently extensive 
to prove this effect [[2 3 . The parameters of equation (12) are 
evaluated at a film temperature defined by equation (1), The use of a 
film temperature allows for variations of all physical properties and 
resulted in the best correlation by an empirical equation. 

For equal capacitance laminar runs, the inner wall temperature 
correction for thermocouple location was made using Martinelli' s ([ibj 
turbulent temperature distribution rather than the laminar temperature 
distribution. This was done because it gave the best correlation of 
these data /ith respect to constant wall temperature data and the turbu- 
lent equal capacitance data. 

Figures 11, 12, 13, and 16 are dimensionless plots of annulus wall 
temperatures versus annulus length, and show that for the flow conditions 
of this investigation the temperature distributions are definitely not 
linear. However, it is believed that the arithmetic mean of the fluid 
bulk temperature is a close approximation to the average fluid tempera- 
ture for equal capacitance runs. The arithmetic mean bulk temperature 
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for the constant wall temperature runs was found to vary as much as 
25% from the average bulk temperature as found by applying Simpson's 
Rule. It is recommended that the average bulk temperature be used to 
evaluate the heat transfer as well as the film temperature for evalua- 
ting fluid properties. The average bulk temperature for the constant 
wall temperature runs was found to be 

=• . 4»73 4 tbct ^ ^2) 

and is recommended for use where knowledge of the bulk temperature dis- 
tribution is unknown. 

All of the laminar runs were also correlated using the Chen^ Hawkins 
and Solberg equation (10) but basing the fluid properties on a film tem- 
perature and using the above average bulk temperature, equation (13). 

These data plot an average of 15% higher than predicted by equation (10). 

A plot is shown in Figure 9. 

A final correlation was made for all runs based on the Colburn an- 
alogy [19] of 

(14) 

with the fluid properties evaluated at the film temperature, equation 
(1). These results are shown in Figure 17. While It was expected that 
the Colburn analogy would produce an acceptable curve for the fixed 
geometry of the system and counter flow fluid arrangement, it was indeed 
startling to find the small scatter of the data and the straight line 
behavior exhibited. It can only be postulated from the correlation that 
for the Reynolds number range covered and the small temperature differences 
which were employed, natural convection effects for all of the runs must 
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have been of the same order of magnitude. Any attempt to introduce 
the Grashof number into the correlation merely scattered the data^ thus 
indicating the free convection effects were negligible. The extrapola- 
ted conclusion can be only that for a certain range of diameter ratios 
and equivalent diameters, even for laminar flow conditions, the Grashof 
number is not a pertinent parameter. 
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VII, Conclusions and Recommendations. 



The friction data Indicate that the proper Reynolds number to 
employ Is c/« , where de Is the equivalent diameter of the annulus. 

The friction data In the laminar region check the theoretical relation 
of Lamb C 7 ] and this should be employed for annuli In this region of flow. 

The heat transfer data, both laminar and turbulent, were correlated 
by equation (12) based on the average bulk temperature and the fluid pro- 
pettles evaluated at a film temperature defined by equation (1). This 
equation Is presented as being valid only for use In Investigating this 
particular geometry annulus and for Identical conditions of flow. It 

■a 

is also suggested that equation (12) Is valid for both heating and cool- 
ing. 

There is to date no single empirical correlation which can account 
for all the possible variationsof flow conditions and annuli geometry. 

If knowledge of the heat transfer within + 57 , is needed, then a particu- 
lar correlation must be found in the literature that applies to the 
specific type of flow pattern and annuli geometry under investigation. 

It would be interesting to know exactly how natural convection 
effects the heat transfer in vertical annuli. Tills could be found by 
reversing the direction of flow In the annulus and tube and analyzing 
Che effect. It has to be expected that for large Reynolds numbers and 
small Grashof numbers, the influe,nce of free convection on the heat trans- 
fer can be neglected. At what point this occurs is yet to be determined. 
The Grashof number varied as much as 200% over the range of Reynolds 
numbers investigated. The Grashof number enters into equation (10) to 
the 0.05 power only and as such varied less than 2% over the range of 
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Reynolds numbers investigated. However, if the Grashof number enters 
into the correlation to the 0.065 power as in equation (12), it will 
produce a 12% variation over the range of Reynolds numbers investigated. 
This agrees with Eckert and Diaguila's CtG prediction of a 10% varia- 
tion in the heat transfer due to natural convection effects. An account 
must be made as to how the buoyancy forces act in conjunction with the 
external forces. A variation of 10% in the Nusselt number would result 
from ignoring this fact. 
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Table I. Apparatus Dimensions 



Inner Tube: 1 inch I. D. x 20 foot Copper Tube 

Effective tube length - 142 inches - 11.82 feet 
Diameter (inside) - 0.995 inches - 0.0825 feet 
Diameter (outside) - 1.125 inches - 0.0948 feet 
Wall Thickness - 0.065 inches 

Cross-Sectional Area (inside) - 0.00498 square feet 
Circumference (outside) - 0.29437 ft. 

Effective Heat Transfer Surface (outside) - 3.48 square feet 
Jacket Tube: 

Diameter (inside)- 1. 242’ inches - 0.1035 feet 
Diameter (outside) 1,3^2 inches - 0.JH3 feet 
Cross Sectional Area (inside) - feet 

Circumference (inside)- 4.46 inches - 0.372 feet 
Annulus : 

Equivalent Diameter; 0.117 inches - 0.00976 feet 

Cross-sectional Area - 0.00151 feet 

Effective friction length - 100 inches - 8.33 feet 
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TABLE II 



Isothermal Friction Data 
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196 
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0,00932 


7960 
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TABLE III 

Constant Wall Temperature 
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